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Abstract 



The drag reduction properties of a turbulent channel flow modified by spanwise 
sinusoidal oscillations of the walls are investigated by direct numerical simulations. 
The work is based on the linear relation between the drag reduction and the param- 
eter S, function of the maximum wall velocity and the period of the oscillation. This 
quantity, first determined by Jl-Choi et al. (2002) and later studied by Quadrio and 
Ricco (2004), has been found through physical arguments pertaining to the action 
of the oscillating Stokes layer on the near- wall turbulence dynamics. The predictive 
potential of the scaling parameter is exploited to gain insight into the drag-reducing 
effects of the oscillating wall technique. The period of oscillation which guarantees 
the maximum drag reduction for a given maximum wall displacement is studied 
for the first time. The issue of the minimum intensity of wall forcing required to 
produce a non-zero drag reduction effect and the dependence of the drag reduction 
on the Reynolds number are also addressed. The drag reduction data available in 
the literature are compared with the prediction given by the scaling parameter, thus 
attaining a comprehensive view of the state of the art. 

Key words: Turbulent channel flow, turbulent drag reduction, spanwise wall 
oscillation, direct numerical simulation 



1 Introduction 



This paper presents a numerical investigation of a turbulent channel flow 
with sinusoidal spanwise oscillations of the walls. The flow over the oscil- 
lating walls results from the combination of two simpler flows, i.e. a canon- 
ical turbulent channel flow in the streamwise direction and an oscillating 
boundary-layer-type flow in the transversal (spanwise) direction. The most 
relevant characteristic of this modified turbulent flow is the time- sustained 



reduct ion of the streamwise wall-shear stress, first pointed out by I Jung et al 



( 119921 ). The drag reduction effect is caused by the weaken ing of the rele 



vant t urbulence-producing events in the vicinity of th e wall (jAkhavan et al 



1993 : iDhanak and S i. 199 91: IChoi and Clavtonl . l200ll ; jPi Cicca et al 



2002 



Iuso et all . 120031 : iRiccol . 12004 ; Ku and Huand . 120051 : IZhou and Balll . |2006|) 
but th e precise details are still poorly understood. I t has also been estab 



lished (IBaron and Quadrid 119961 ; iQuadrio and Riccol . 120041 ) that a net en- 



ergetic saving of the order of 10% (determined by taking into account the 
power spent to move the walls against the viscous resistance of the fluid) 
can be obtained by carefully tuning the parameters of the oscillation. An- 
other property of this flow that will be exploited in the following is that, once 
the flow field is averaged along the streamwise and spanwise homogeneous 
directions, the spanwise velocity pr ofile agrees with the laminar solution of 



the so-called second Stokes prob lem (jQuadrio and Sibillal . l2000l ; iJI-Choi et al. 



20021 ; IQuadrio and Riccol . 120031 ) , hence uncoupling from the complex dynam- 



ics of the all-encompassing turbulence. In the spirit of the wall-oscillation 



ing the turbulence by spanwise travelling waves (Du and Karniadakis. 2000; 


Du et al.. 2002: 


Karniadakis and Choi. 


2003; 


Zhao et al.. 2004: Itoh et al.. 


2006 


Yoon et al. 


2006). bv spanwise oscillating Lorentz forces (Berger et al.. 


2000 


Pane; and Choi. 2004; Breuer et al. 


2004: 


Lee and Sung. 2005). bv span- 



1/ I >. f tl 1 1/ 1 ^i-^ ^^M^^H-W 1 1/ / U ± 

wise oscillating suction and blowing (ISegawa et all 120051). and by stead y 



streamwise oscillations of the spanwise wall velocity (jQuadrio et all 120071 ). 
The wall oscillation has also been shown to be effective in reducing the growth 
rate of the most unstab le Gortler vortex developing on a concave surface 
(iGalionis and Hall l2005h . 



An important step towards practical applications of the oscillating-wall tech- 
nique is the recent finding of a scaling parameter S that is suggested to be 
related to the amount of drag reduction. This parameter depends on the quan- 
tities defining the sinusoidal oscillation, namely the period of oscillation T and 
the maximum wall velocity W m . (A third parameter describing the oscillation 
is the peak-to-peak wall displacement D m , which is D m = W m T/ii for a si- 



nusoidal waveform.) Jl-Choi et al.l ( 20021 ) have correlated their drag reduction 
direct numerical simulation (DNS) data with S, which was found from phys- 
ical arguments pertaining to the interaction between the spanwise laminar 
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Stokes layer and the n ear- wall turbulence. A leas t-squares fit yielded a power 



law expression for S. iQuadrio and Riccol (|2004f) f den o ted b y QR hereafter 



have recently improved the analysis by iJI-Choi et al.l (120021 ) on the basis of 
a DNS dataset of a turbulent channel flow modified by the motion of both 
walls, for a Reynolds number Re T = 200 defined by the friction velocity of the 
natural turbulent flow and half the distance between the channel walls. The 
accuracy of their dataset brought to light a linear relation between the drag 
reduction and S, which was suggested to hold as long as T remains smaller 
than a typical integral time scale of the wall turbulence. 

The main objective of the present work is to investigate further the possi- 
bility of estimating the drag reduction through the parameter S. In QR the 
correlation analysis between S and the drag reduction was carried out only 
with their DNS data. In the present analysis all the drag reduction datasets 
available in the literature are employed to compute a new, and more general, 
least-squares fit. A useful by-product of this procedure is that each dataset 
(both from experiments and DNS) can be discussed in terms of its difference 
from the fit. Since the literature data show a considerable scatter, comparing 
the actual measurements of drag reduction with the predictions via S permits 
an indirect assessment of their accuracy. 

The linear fit allows us to identify and study easily two distinct periods of wall 
oscillation that are optimal in terms of drag reduction. Although in DNS stud- 
ies the oscillation parameters can be chosen freely, T optjW , the optimal period 
of oscillation for fixed maximum wall velocity, has usually been the quantity 
of interest. On the other hand, an experimental campaign is likely to reveal 
T op t,D, the optimal period at fixed maximum wall displacement. Indeed, in a 
laboratory the wall motion is usually produced by a crank-slider mechanism, 
which allows varying T for fixed D m . QR however noted that T opt ^ had never 
been observed by experimentalists, probably because the flow at these high 
frequencies of oscillation is difficult to test. The difference between T opt jj and 
T op t,w was left unnoticed until QR pointed out that T opt ,D should be smaller 
than T optt w, which is known to be constant at T^ pt w rs 125 (when scaled by 
the friction velocity of the undisturbed flow, Karniadakis and Choi, 2003). 
QR's analysis also showed that T opt D should be a function of D m itself. Al- 
though their limited data appeared to validate these unconfirmed conjectures, 
it is an important aim of the present study to verify the existence of T optj D 
and to assess its properties. 

The analytical expression for S yields additional information. As a further 
aim, we shall use S to seek the smallest values of the oscillation parameters 
which guarantee a non-zero drag reduction effect. These oscillating conditions 
are referred to as "minimal" throughout the paper. The same analysis will also 
allow us to determine the oscillatory conditions needed to attain a specified 
amount of drag reduction and of net energy saving. 
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We shall also address the important issue of the effect of Reynolds number on 
the drag reduction. This point is still open to discussion, as the conclusions 
of previous studies do not fully agree. The ex perimental analysis of a turbu- 



lent boundary layer by iRicco and Wul (120041 ) has shown that increasing the 



Reynolds number based on the momentum thickness and free-stream velocity 
from Rep =500 to 1400 ( with ~ 240) had no influence on the drag reduc- 



tion. I Choi and Graham! (Il998l )'s experimental results in a pipe flow modified 
by circular wall oscillations have also indicated that changing the Reynolds 
number based on the bulk velocity and the pipe diameter from i?e=23,300 to 
36,300 gives variations of drag reduction which are of the order of the uncer- 



tainty range. The DNS of a turbulent channel flow conducted by lJI-Choi et al. 



( 120021 ) has instead revealed that the reduction in wall friction can be halved 
by increasing the Reynolds number from Re T = 100 to Re T = 400. We shall 
carry out a few DNS to ascertain whether or not increasing Re T affects the 
drag reduction. 

The paper is organized as follows. In Section 2, the numerical procedure and 
the discretization parameters are described. The main steps of the analysis by 
QR are briefly recalled in Section 3, where the parameter S and the function 
Topt,D = To P t,D(D m ) are defined. Section 4 discusses the quality of the analytical 
prediction of T opt! D on the basis of the DNS results. The minimal wall forcing 
conditions and the dependence of the drag reduction on the oscillatory param- 
eters are presented in Section 5. Section 6 describes the effect of the Reynolds 
number on the drag reduction. Section 7 presents the drag reduction data 
available in the literature, providing an evaluative overview on the state of 
the art in terms of a comparison between the amounts of drag reduction and 
the estimates based on S. Section 8 is devoted to a summary. 



2 Numerical procedure 



We have studied the turbulent flow in a channel with moving walls through 
direct numerical simulations of the incompressible Navier-Stokes equations. 
The walls move in phase along the spanwise direction with velocity 

W(t;W m ,T) = W m sm (^-Y 



Figure 1 displays a sketch of the computational domain. 



Our pseudo-spectral solver is described in iLuchini and Quadrid (120061 ): it is 
based on Fourier expansions in the homogeneous directions and on fourth- 
order accurate, compact finite-difference schemes for the discretization of the 
differential operators in the wall-normal direction. Aliasing errors in the com- 
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putation of the non-linear terms are eliminated by expanding the flow variables 
into a (at least) 3/2 larger number of modes for each homogeneous direction 
before transforming from the Fourier to the physical space. The temporal 
integration is carried out with a partially-implicit procedure: a third-order, 
low-storage Runge-Kutta method for the convective terms, and a second-order 
Crank-Nicolson scheme for the viscous terms. The mixed spatial discretization 



is adv antageous from the viewpoint of parallel computing (ILuchini and Quadrio 



20061 ). and allows employing computing machines connected by standard net- 



working hardware to achieve a large computational throughput. 

The simulations described in this paper have been carried out on a computing 
system available in dedicated mode at the University of Salerno, made by 64 
machines, each with two Opteron CPUs. We have performed the calculations 
at two values of the Reynolds number: Re T = 200, 400 based on h, half 
the distance between the channel walls, and on u T , the friction velocity of 
the uncontrolled case. The computational parameters for Re T = 200 have 
been chosen to replicate those in QR, if exception is made for a slightly finer 
wall-normal discretization. We have employed a computational box with the 
following dimensions: L y = 2h, L x = 21h, L z = 4.2h in the wall-normal, 
streamwise and spanwise directions, respectively. We have used 160 grid points 
in the wall-normal direction, and 321 and 129 Fourier modes in the streamwise 
and spanwise directions. For the Re T = 400 case, we have employed 256 grid 
points in the wall-normal direction, and 385 and 321 Fourier modes in the 
streamwise and spanwise directions. Each simulation has been run for 40,000 
time steps with At + = 0.2. Throughout the paper, the + superscript indicates 
quantities scaled by the inner units of the unperturbed turbulent flow. A few 
calculations with very high W m have required an even smaller (up to one half) 
time step, due to stability constraints. A single case at Re T = 200 has taken 11 
seconds for one time step and 5 days for the full simulation with 8 machines. 
The wall-clock time for the whole set of simulations has been slightly more 
than two weeks with the full use of the computational power. 

For the Re T = 200 case, we have carried out one simulation for the canonical 
flow and 28 simulations for different (T, W m ) pairs. The skin-friction coefficient 
for the fixed-wall case is Cf = 2r x /(pU b 2 ) = 7.94 x 10~ 3 , where r x is the time- 
and space-averaged streamwise wall-shear stress, p is the fluid density and Ub 
is the bulk velocity. This value essentially coincides wit h the value estimated 



by the following empirical formula given in iPopd (120001 ) at page 279: 



Cf = 0.0336/ter - 273 . 



Three values of D m have been considered, = 100, 200, 300. At a given D m , 
we have explored the existence of T optt D by varying T along the hyperbola 
W m = D m 7i/T in the (T, W m ) space. The amount of drag reduction has been 
determined by a procedure (QR) which involves discarding the initial temporal 
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transient, averaging over time and along the homogeneous directions, and 
considering both walls to enlarge the statistical sample. For the Re T = 400 
case, we have carried out one simulation for the canonical flow (Cf = 6.495 x 
10~ 3 , only 1% smaller than the value given by (1)) and three simulations with 
W+=12 and T+=30, 125 and 200. 



3 Laminar analysis 



The analysis is based on the close agreement between the space-averaged 
(along the homogeneous directions) spanwise flow and the l aminar solution for 
the second Stokes problem , defined by wt (y + , t + ; T + ) (IQuadrio and Sibillal . 
2000l : Ijl-Choi et all . 120021 ; buadrio and Riccol . 120031 ). At relatively high fre- 
quencies, the space-averaged spanwise flow uncouples from the space-averaged 
streamwise flow, although the spanwise turbulent fluctu ations are significantly 
altered by the wall motion ( Quadrio and Riccol . 2003 ). We can explain this 
agreement as follows. The space-averaged spanwise momentum equation is 



dw + 
~di+ 



d 2 w 
dy 



+2 



— R v 



t+ 



where R vw = dv'w' + /dy + \ t +. The barred quantities are averaged along x and 
z, w is the turbulent spanwise velocity, and v' and w' are the wall-normal and 
spanwise velocity components fluctuating about their corresponding space- 
averaged quantities. The equation describing the laminar Stokes flow is ob- 
tained by replacing W with w s , and by setting R vw = 0. Fig. 2 shows R vw at 
different oscillation phases from the start-up of the wall motion for = 18 
and T + = 125. It initially increases, reaches a maximum when t ~ T/2, and 
then decreases as the transient elapses after a few oscillation periods. As R vw 
is negligible once the new modified turbulent flow is established, it follows that 
w + « as t + — > 00. The fact that w + does not coincide with during the 
initial transient is reflected in the time history of the spac e-averaged spanwise 
wall-s hear stress r z . It is observed in Fig. 8 at page 12 in Quadrio and Riccol 



(120031 ) that this quantity deviates slightly from the laminar value at t ~ T/2, 



when R mi , reaches its maximum value. 



The sc aling parameter S is now introduced. It has been expressed by lJI-Choi et al 



( 120021 ) as the product of a+ , the maximum acceleration of the Stokes layer at 
a specified wall-normal location, and a distance from the wall at which the 
oscillating layer directly affects the turbulent flow. 
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At the ro ot of this proced ure lies the laminar solution for the second Stokes 
problem (jBatchelorl . 119671 ) : 



2tt 



W£exp ( -y + Jn/T+ ) sin ( ^-t + - y + J^- ) .(2) 



The maximum spanwise acceleration a+ at a wall-normal distance y + = is 
obtained by differentiating (2) with respect to t + : 

«i=^«p(-<:^). 



It is further required that the maximum spanwise velocity of the Stokes layer 
at y + = be larger than a typical value of the spanwise velocity fluctua- 
tions. The wall forcing must be intense enough for the Stokes layer to influence 
the turbulent fluctuations and disrupt the near-wall turbulence-producing cy- 
cle. By imposing the above condition and by using (2), it follows 




The quantities a+ and are united into the expression for S: 



and have been determined by maximizing the correlation coefficient 
Cs between the drag reduction data and S. In QR, C$ = 0.99, = 6.3 and 
Wfa = 1.2 (which remarkably agrees with t he maximum r.m.s . of the turbulent 



spanwise velocity fluctuations, w+ ms wl.l (IKim et all 119871 )) and 



DR(%) = S 1 S + S 2 , 



(4) 



where DR(%) is the drag reduction, 5i = 130.6 and 5*2=-2.7. 

An expression for T^ ptD , the period of oscillation which guarantees the maxi- 
mum drag reduction for fixed D m , is found by setting: 



dDR 
dT+ 



Si 



D2 



OS 
df+ 



0. 



IK 



d 2 DR 
<9T+ 2 



Si 



d 2 S 
8T+ 2 



< 0. 
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The second condition has been verified graphically via (3). After eliminating 
W+ = ttD+/T+ in (3), it follows: 




1 In 



7T 



D: 



, L o V t,D™ th, 



2. 



(5) 



Differently from Tt^, which does not depend on W+, D increases mono- 
tonically with .D+ and is smaller than T^ t w . The latter is obtained by dSj_ dT 
0, i.e. T^ vt w = Ti (£+ ) 2 ~ 12 5 , which agrees with the v alues in lJung et al. 



(119921 ) . iDhanak and Sil (I1999h . buadrio and Sibillal (J200rJ) and QR. 



QR have also shown that S scales linearly with the drag reduction only for 
T + < 150. Fortunately, it occurs that T^ t D ,T^ ptw <150, so that the pre- 
diction of the optimal periods based on the linear relation between the drag 
reduction and S is valid. These quantities are not well correlated when the 
oscillation uncouples from the near-wall turbulence dynamics, namely when 
half the period of oscillation is larger than a typical pseudo-Lagrangian time 
scale r epresenting the survival time of the longest-lived structures (~ 60 time 
units) (jQuadrio and Luchinil . 120031 ) . When T is large, the near- wall structures 
have enough time to develop their inner dynamics between successive sweeps 
of the Stokes layer. In this limit, the flow adapts to a new quasi-steady three- 
dimensional condition, where time can be treated as a parameter and the 
drag-reducing effect of the Stokes layer is lost. 



4 Drag reduction scaling and optimum period at fixed D r 



Fig. 3 shows the drag reduction DNS dataset versus S for Re T = 200 produced 
for the present analysis, together with the data by QR for T + < 150. The 
oscillating conditions are largely different, but the data collapse well on the 
straight line. The linear regression (4) has been recomputed by grouping our 
DNS data with most of the data available in the literatur e, which are analysed 
in Section 7. The only discarded datasets were the ones by lJI-Choi et al.l (120021 ) 
at Re T = 100 becau se during the oscillatio n the Reynolds number was too 
low and the ones by IDhanak and Sil (Il999l ) , whose analysis was based on a 
simplified model and not on the Navier-Stokes equations. The new correlation 
parameters are C s = 0.92, 5i=135.11, S 2 =-0.85, = 6.2 and W t + h = 1.7, 
which are essentially unchanged from the analysis in QR for Re T = 200. These 
new values have been used for the present data analysis. 



The results of the present simulations, designed to identify T op t t D, are com- 
pared in Fig. 4 with the prediction obtained by (3) and (4). Good agreement 
occurs except for three data points at low T (open symbols in Fig. 4 and black 
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squares in Fig. 3). For these points, the predicted values are lower than the 
actual DNS data. We have not been able to explain this behaviour at very 
small T, which corresponds - being the displacement fixed - to very large W m 
(W+ >40). We have first verified that the space-averaged spanwise velocity 
profile at various phases still agrees with the laminar solution, from which (3) 
is determined. As a further check, the Reynolds number Res for the spanwise 
oscillating flow, based on W m and on the Stokes layer thickness 5 S = yj vT/n 
(where v is the kinematic viscosity of the fluid) is compared with the criti- 
cal Reynolds numbers for stability or transition of the Stokes flow, although 
this analysis is obviously not rigorous in that we consider a space-averaged 
profile and not a purely laminar flow. Despite the wall velocity being as high 



as twice the centreline velocity, Res = W+yT + /n is not high, since T is 
low. For our data, Res = 8 — 1 40, which is lower than « 1400 found by 



Blennerhassett and Bassoml (120061 ) as the lowest critical Reynolds number for 



an oscillatory boundary l ayer b etween two parallel plates, and ~ 550 reported 
by IVittori and Verziccol (Il998l ) as their transition threshold. 



However, these discrepancies do not limit the possibility of employing S to 
study successfully the flow because all the relevant drag reduction features 
occurs for smaller W£. Indeed, (i) T optiD (D m ) approaches a constant value 
for >25 (see Fig. 1 at page 259 in QR), (ii) the explored range of wall 
velocities has been <16 in the previous experimental works and W£ <18 
in the previous numerical works (except for a few cases at W+=27 in QR), 
(iii) QR have found that the drag reduction does not change significantly with 
W m for >20, which greatly simplifies the analysis and it implies that an 
analysis through S is not needed at these high wall velocities, (iv) it will be 
shown in Section 5 that the net energy saving may be positive only in the 
range <7. 



The fact that the prediction for T opt)D is good over most of the explored range 
is further confirmed by Fig. 5, which shows that the optima computed by the 
DNS data compare satisfactorily with T opt D = T opt D (D m ) obtained by (5). 



5 Minimal oscillating conditions and estimates of DR{%) and P ne t(%) 
as functions of W+, D^, T + 



The idea that a finite intensity of the forcing is needed to affect the turbulent 
friction is contained in the definition of S, where a threshold velocity Wth is 
introduced. The fact that the regression line in Fig. 3 crosses the abscissa at 
S = S m i n = 0.0063 > implies that the wall must oscillate with a minimal 
velocity W m , min > W th (or with a minimal displacement D m ^ min > W th T/n) 
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to achieve drag reduction. From (3) it follows that 



| ' ' m,min 


> = < 




\d+ . 

^ m,min 







W t + h exp 



Q . l r F-\ 



7T 



exp 



7T 

T+ 



(6) 



These minimal quantities are displayed in the contours plots of Fig. 6 as the 
zero drag reduction curves. The drag reduction for a (T + , W£) pair, computed 
by (3) and (4), is represented in Fig. 6 (top graph) for 30 < T + < 150. W m ^ min 
grows unbounded as T decreases because a stronger wall forcing is needed 
to affect the turbulent flow as the penetration depth 5 S ~ \/T of the Stokes 
layer vanishes. W^ imin becomes approximately constant at higher periods of 
oscillations, say T + > 30, and its value ~ 1.8 is of the order of magnitude of the 
near- wall spanwise velocity fluctuations. It also follows that D^ min ~ D{T + , 
D 1 = W+ min (T + = 150) /tt = 0.57, for 30 < T+ < 150. 



Fig. 6 (bottom graph) shows the drag reduction corresponding to a (T + , D+) 
pair. Dtn,min(^optyv = 125) ~ 70 compares well with the spanwise width of a 
low-sp eed streak, pa 20— Q y/u T from flow visu alizations of the near- wall turbu- 
lence ( Hirata and Kasagl 1979 ; Ricco, 2004 ). This result confirms that these 
structures need to be suffi ciently swept laterally to achieve drag reduction 
( Baron and Quadriol . 1996 ). The minimal conditions are in good agreement 
with the experimental results by Raskob & Sanderson (personal communica- 
tioJD), who first showed that a displacement w 80 is needed for drag 
reduction when 80 < T + < 1000. The existence of the minimal conditions is 



also confirmed by the analysis of iQuadrio and Riccol (120031 ). They have found 
that the wall-shear stress is not affected at the beginning of the oscillation, 
when the wall velocity is small. Figs. 2 and 3 in their paper show that the 
space- averaged r x changes by less than 1% when W + < 2. This value com- 
pares well with W^ min and it remains essentially unchanged as T + varies for 
50 < T + < 200, similarly to the behaviour of W^ min in the same T + range 
(see Fig. 6, top graph). 



An an alogous behaviour emerges from t he experimental stu d y bvlMao and Hanrattv 
( 1986h and by the numerical works by Ismael and Cotton ( 1996 ) and Cotton 
(l2007h . who studied the variation of wall-shear stress in turbulent pipe flows 



subjected to streamwise oscillations of the pressure gradient. They found that 
for small forcing amplitudes the flow responds linearly and the mean value of 
the velocity gradient at the wall is unaffected. 



B. Raskob and R. Sanderson presented these results at the APS Division for Fluids 
Dynamics Meeting in San Francisco, California in 1997 with the title Turbulent drag 
reduction due to an oscillating cross-flow. 
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The expressions (6) contain the parameter S m i n . Table 1 shows that this quan- 
tity varies substantially when and £^ change by a small amount. However, 
the minimal conditions remain almost unvaried, as revealed by the values of 
W^ min at T+ = 125 in Table 2. 

We now turn to the prediction of the net power saving P ne t, defined as the 
difference between the power saved thanks to the wall motion, i.e. the drag 
reduction when [/& is constant, and the power P sp spent to move the walls 
against the viscous stresses. These quantities are expressed as percentage of 
the power spent to drive the fluid through the unmanipulated channel, 2t^U^~ , 
where = 15.88 at Re T = 200. P sp over an interval tf — U is 




where I and u denote the lower and upper walls. Inasmuch as r z (t) is well 
predicted by the laminar Stokes solution (see Section 3), (7) becomes 



P SP (%) = 



100 



T+ 

I 



w~ 



dwf 



dt 



2U+ 



IT 



T+ 2U+Re s 



,(8) 



P net (T + ,W+) and P net (T + ,D^ n ) are displayed in Fig. 7. Steep changes occur 
when P net < 0, whereas regions with P net > present more gradual variations. 
When compared with DNS results by QR, the period T + 60 at = 4.5 
needed for P net > is well captured. For T+ = 150, P net > for 2 < < 7, 
while in QR this interval is 1.5 < < 9. 



By setting 

dP ne t 
dW+ 



o, 



T+ 



d 2 p 



net 



dW+ 2 



<0, 



T+ 



the maximum wall velocity and displacement which give P ne t,rnax, the maxi- 
mum P net at fixed T + , are found: 



' " net,max 
^net,max 




>exp l-f 



7T 

T+ 



where W x = = 6.55 for Re T = 200. It follows that 

Pnet,max{T + ) = S 1] f^e W ("C^) (A " % \ff*) + ft, 
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where Pi = 2\n(W 1 /W t t) -1 = 1.7 for Re T = 200. The location of the overall 
maximum P net is well estimated at T + = 150 and = 4.2, while its value 
is slightly underpredicted: 5.6% instead of 7.3% (QR). 



6 Effect of Reynolds number 



Our calculations show that slightly lower amounts of drag reduction are ob- 
tained by doubling the value of the Reynolds number from Re T = 200 to 
Re T = 400 for W+ = 12 and T+ = 30, 125 and 200, as shown in Fig. 8. 
The wall-shear stress reduction decreases from 21.7% to 20.3% at T + = 30 
(6.6% change), from 32.5% to 28.1% at T+ = 125 (13.4% change), and from 
27.2% to 22% at T+ = 200 (19.2% change). These variations thus increase 
with T for fix e d Wm . This result is in broad agreement with the analysis by 
Jl-Choi et al.l ( 20021 ). although a smaller computational domain and a lower 



spatial resolution were employed in their analysis. Their values of drag reduc- 
tion at Re T = 200 are lower than ours and their drag reduction variations 
with Re T at = 10 are more significant, namely 11.2% at T + = 50, 25.1% 
at T+ = 100, 29.9% at T+ = 150, and 29.3% at T+ = 200. 

Previous experimental studies indicate that the drag reduction does not vary 
with the Reynolds n umber for (low) values of this parameter. For example, 
Ricco and Wul (120041 ) have found that in a turbulent boundary layer the drag 
reduction does not change as the Reynolds number based on the momen- 
tum thickness 9 and the free-stream velocity increases from Peg =500 
to 1400 (which co r respo nd^] to Re T =257 and ite T =633), when L>+ « 240. 
Choi and Graham fjl998h have shown that in a pipe flow modified by circular 
wall oscillations the wall-shear stress reduction remains within the range of the 
experimental uncertainty when the Reynolds number Re based on the bulk 
velocity and the pipe diameter increases from i?e=23,300 to 36,300 (which 
correspond" 3 ! to Re T = 650 and Re r = 962). The uncertainty range on the per- 



2 The Reynolds number Re$ for a free-stream boundary layer can be converted to 
Re T = 5u r /u (where the boundary layer thickness 5 is the wall-normal distance 
where the mean velocity is 0.99f/oo) for easines s of comparison wi t h the channel 
and pipe flow data. By assuming that 5 = 1W (Bogard and Thole] . 19981 ) and by 



using Cf = 2(u T /U OD ) 2 = 0.025i?e^ a25 (|Kavs and Crawford! . Il993l ). one arrives at: 



Re T = 1.118-Re; 



0.875 



Re e < 3000. 



(9) 



3 The Reynolds number Re for pipe flow can be co nverted to R e T based on u T and 
pipe radius by using the following formula given in Pope ( 2000l ) at pages 292-293: 



Re = iV2Re T [21og 10 (iV2Re T ) - 0.8 



(10) 
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cent drag reduction was about 10% for IRicco and Wul (120041) (see T able 1 at 
page 45 in their paper), and about 20% for Choi and Graham ( 19981 ) (see Fig. 
2 at page 8 in their paper). It was higher in the second case probably because 
the wall-shear stress has been determined via measurements of the mass-flow 



rate and pressure drop, while IRicco and Wul (120041 ) have measured directly 
the mean streamwise velocity in the viscous sublayer. It is thus possible that 
the drag reduction variations have not been detected because of the high ex- 
perimental uncertainty. This conjecture is further s upported by the fact that 
our numerical uncertaint y is lower than the ones by IRicco and Wul (120041 ) and 
Choi and Grahaml (119981 ) on account of the difficulty of such experiments and 
the high accuracy of our computations (refer to section 2 .4 at page 256 in QR) . 
We also note that the majority of the experiments by IRicco and Wul (120041 ) 
have been conducted at relatively low T + (T + pa 50, D+ pa 240, W+ pa 15 - 
see Fig. 13 in their paper), where the change of the drag reduction with the 
Reynolds number is weak, as shown in Fig. 8. 



Our results might also explain why numerical works have revealed higher drag 
reduction values than the experimental ones (see Fig. 14 at page 51 in Ricco 
and Wu, 2004), the latter generally conducted at higher Reynolds numbers. 
The effect of Reynolds number however warrants further investigation as we 
have only explored a portion of the (T, W m ) space. It would be of interest to 
verify the existence of a maximum Reynolds number above which the wall- 
shear stress is unchanged, and to investigate how it varies with T and W m . 

The effect of the Reynolds number on the maximum net energy saving can 
be estimated by considering separately the changes on the power spent P sp 
to move the walls and on the power saved through the wall motion, i.e. the 
drag reduction. As shown by (8), P sp (%) varies with the Reynolds number 
only because of a change of Uj}~, which can be expressed as follows: 



u: 



Re 
2Re T 



7.715Re 



0.136 



where we have used Re = 15 A3Re]. given by iPopd (120001 ) at page 279. (U^ 
at Re T = 200,400 obtained by our calculations is less than 1% different from 
the one given by (11)). It follows that 



P sp (%) = 6.481(W+)V ^e; ' 136 . 



(12) 



The power spent for = 4.5 and T + = 125, namely for the oscillating 
conditions at which the maximum net energy saving occurs (see Section 5), 
decreases from P sp =10.1% at Re T = 200 to P sp =9.2% at Re T = 400. By as- 
suming that the drag reduction for these oscillation conditions decreases by 
the same amount of the drag reduction at = 12 and T + = 125, the max- 
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imum net energy saving is estimated to decrease slightly from P n et,max=7 -1% 
for Re T = 200 to 5.7% for Re T = 400. We have assumed that P ne t,max does 
not change location in the (T + , W*) space, which remains to be verified. Such 
change would only be due to a shift of the drag reduction peak, i.e. T+ t w , and 
not to P sp , which only changes in magnitude when Re T varies, as indicated 
by (12). The data in Fig. 8 suggest that T^ vt w should not vary much as Re T 
increases. 



7 Analysis of available datasets 



Many drag reduction data, either from DNS or experiments, have been re- 
ported in the literature for the oscillating wall technique and they are all af- 
fected by various sources of error. Accuracy concerns for DNS-based datasets 
arise mostly from issues related to the spatio-temporal discretization and the 
computational procedures. In laboratory experiments hot-wire anemometry 
has often been employed and the measurements may have been biased by 
errors due to the highly three-dimensional flow field. A comparison among 
all the drag reduction data have been carried out by QR in the space of the 
parameters (T, W m ), but the analysis did not convey the desired information 
about the accuracy because of the high scatter of the data. We attempt here to 
attain a clearer picture of the state of the art by comparing the S-based pre- 
diction of drag reduction with the wall-shear stress measured in the laboratory 
or computed via DNS. 



Fig. 9 show s a first group o f DNS results. The three dataset s bvlBaron and Quadrio 



1996) and lJI-Choi et aJj (120021 ) for a channel flow and by lQuadrio and Sibilla 
(l2000l ) for a pipe flow are all computed at Re T = 200 and all broadly agree 



with one another. Indeed, the numerical accuracy is similar in the three cases. 
The disparity between the data and the prediction is of the same order of their 
accuracy (which is lower than the one of the present analysis) i n terms of the 



percen t drag reduction. The best correlated dataset is the one by lJI-Choi et al 
fl2002h . 



Jl-Choi et al.l (120021 ) 7 s data in Fig. 10 show higher drag reductions for lower 
Re T . Almost all the data points at Re T = 150 for the pipe flow fall between 
the channel flow data at Re T = 100 and Re T = 200, which may indicate 
a negligible influence of the flow geometry. All the data with T + = 200 are 
below the linear trend, confirming that the S'-scaling occurs only for T + < 150 
(QR). 



In Fig. 11, t he DNS studie s by I Jung et al.l (119921 ) for a channel flow at Re T = 
200 and by iNikitinl fj2000h for pipe flow at Re T = 147 (Re = 4000, see for- 
mula (10)) employed the lowest spatial resolution and the shortest integration 
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time. The comparis on with the ^ -estimates is not satisfactory for I Jung et al. 



( 199211's da t a, wh ile iNikitinl (|2000|)'s points fall slightly below the straight line. 
Jung et all (119921 ) show a 10% increase of wall fri ction for T + 500, m 13, 



which has not been confirmed by other works. iDhanak and Sil (119991 ) 's data 
do not follow the expected trend probably because their analysis is based on a 
simplified model and not on the full Navier-Stokes equations. However, their 
analysis was not aimed at an accurate calculation of the wall-shear stress, but 
at improving the physi cal understand i ng. D espite the model limitations, they 
predicted T+ tw » 90. iMivake et ail (ll997l Vs point at Re T = 150 is slightly 



lower than expected. They might have computed a higher drag reduction (and 
thus attain a better agreement with our prediction), had they continued the 
simulation for a longer time (see Fig. 1 at page 203 of their paper, where the 
wall friction appears likely to decrease further). 



Fig. 12 shows that the best correlated experimenta l dataset is the one b y 
Skandaji fll997f ) ( some of these data are published in (ILaadhari et all 119941 )). 
for a free-stream boundary layer at Reg = 770 (Re T = 375 when formula (9) 
is used), Reg = 980 (Re T = 463) and Reg = 1600 (Re T = 711). The points at 
Re T = 711 are at T + > 150 and show the discrepancy due to the wall motion 
uncoupling from the turbulence dynamics. The points with T + < 150 (actually 
T + < 100) show a good correlation with the straight line. This is expected as 
the Re T effect is comparable with the experimental un certainty at low T + . T he 
experimental data for a free-stream boundary layer by lRicco and Wul (l2004h at 
Reg = 500 (Re T = 257), Re e = 950 (Re T = 451) and Re e = 1400 (Re T = 633) 
agree with our prediction. The Reynolds number effect is small as most of the 
data are for 30 < T+ < 70. 



Trujillo et al.l (119971 ) 7 s experimental data for a free-stream boundary layer at 
Reg = 1400 (Re T = 633), presented in Fig. 13, show the correct slope. They 
are however lower than the linear fit, probably because of the bias caused 
by the spanwise component of velocity on the hot-film measurements in the 
proximity of the wall (Trujillo, personal communication; Choi and Clayton 
2001; Ricco and Wu, 2004). The same data have been corrected in iTrujillo 
(119991 ). thus reaching a better agreement with the line (only the corrected data 
with T + < 150 are inc l uded for clarity i n Fig; . 13). The fact that t he data at low 
T by ISkandaiil (Il997h . iRicco and Wul (120041 ) and iTruiillol (Il999h fall near the 
linear regression could be a sign that the two flat-plate geometries, i.e. a free- 
stream boundary layer as in the experiments and a pressure-driven channel 
flow as in our simulations, might have the same or very sim ilar drag reduction 
properties. The pipe flow data by I Choi and Graham! (119981 ) at Re T = 650 and 
Re T = 962 also have the same slope, but are lower than the linear regression, 
proba bly on account of the high values of Re T . The two points by lChoi et al. 
(119981 ) for boundary layer flow at Re T = 549 are higher than expected. 
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8 Summary 



We have presented a study of the drag reduction effects of spanwise wall oscil- 
lations on a turbulent channel flow based on the direct numerical simulation of 
the inco mpressible Nav i er-St okes equations. A sca ling parameter S, fi r st prq - 



posed by lJI-Choi et al.l (120021 ) and later studied by lQuadrio and Riccol (120041 ). 
is further considered and used to improve our understanding of the main prop- 
erties of this drag reduction technique. 

We have discovered that S is an excellent predictive tool for drag reduction for 
W^m — 40 and 30 < T + < 150. This region of the parameters space includes 
most of the published drag reduction data and is of most practical impor- 
tance. The existence of two optimal drag reduction periods, T optt D and T opti w, 
deduced from the expression for S and verified through our DNS database, is 
now clearly established. T opt ,D, the period of oscillation which guarantees the 
maximum drag reduction at fixed peak-to-peak wall displacement D m , pri- 
marily concerns experimentalists, who are forced to vary T with D m constant. 
T op t,Wi the optimum period of oscillation at fixed maximum wall velocity, is 
the quantity typically searched for by numerical investigators. The function 
T op t,D = T p tt D(D m ) has been determined from the expression for S and con- 
firmed via numerical experiments carried out to the purpose. We have also 
established the minimal wall forcing conditions leading to drag reduction and 
obtained plots of the drag reduction and of the net energy saving as functions 
of the oscillation parameters. Despite some discrepancy with the available 
data, such charts might prove useful for the design of a drag-reducing device. 
Further work is necessary to establish how the wall-shear stress changes from 
the unperturbed condition in the limits T + — > 0, oo for /}+ or constant. 

The amount of drag reduction has been found to decrease slightly as the 
Reynolds number varies from Re T = 200 to Re T = 400 when = 12. This 
effect amplifies as the period of oscillation increases. Further study should be 
conducted to investigate how the optimal periods of oscillation change with 
the Reynolds number. We have also estimated that the maximum net energy 
saving decreases slightly with Re T . The existing numerical and experimen- 
tal data have been re-examined through a comparison between the measured 
amounts of drag reduction and their S'-based estimates: differences have been 
discussed on a case- by-case basis in relation to the accuracy of each dataset. 
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Nomenclature 



Symbol Definition 

a m maximum acceleration of Stokes layer, m/s 2 

Cs correlation coefficient in ^-scaling regression analysis, dimensionless 

Cf skin-friction coefficient Cf — 2r x /(pU^) for channel flows 

and Cf — 2T x /(pU^) for boundary layer flows, dimensionless 

D\ coefficient used for prediction of D mj7n i n , dimensionless 

D m peak-to-peak maximum wall displacement, m 

D m ,min minimum value of D m for drag reduction, m 

DR percent drag reduction, (%) 

h half channel height or pipe radius, m 

I superscript indicating lower wall 

l a wall-normal distance at which a m is computed, m 

l m wall-normal distance used in the computation of S, m 

L x streamwise length of computational box, m 

L z spanwise width of computational box, m 

P ne t net power saved to drive fluid through channel 

thanks to wall oscillation, kg m 4 /s 3 

Pnet,max maximum P net at fixed T + at given Re T , kg m 4 /s 3 

P sp power spent to move walls, kg m 4 /s 3 

Rvw Rvw = dv'w' + /dy + , m/s 2 

Re Reynolds number based on Ub and 2h for channel flows 
and on Ub and pipe diameter for pipe flows, dimensionless 

Res Reynolds number for Stokes layer (Reg = W m 5 s /u), dimensionless 

Re e Reynolds number based on U^ and 9, dimensionless 

Re T Reynolds number based on u T and h, dimensionless 

S drag reduction scaling parameter, dimensionless 

Si, 52 constants in equation (4), dimensionless 

S min minimum value of S for drag reduction, dimensionless 

t time, s 

tf time at which computation of wall-shear stress is completed, s 

ti time at which computation of wall-shear stress is initiated, s 

T period of wall oscillation, s 

T opt D optimum period of wall oscillation at fixed D m , s 

Topt,w optimum period of wall oscillation at fixed W m , s 

u superscript indicating upper wall 

u T friction velocity for fixed wall configuration (u T = Jr x /pY m/s 
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Ub bulk velocity in channel and pipe flows, m/s 

Uoo mean free-stream velocity for boundary layer flows, m/s 

v' wall-normal turbulent velocity fluctuating about x — z average, m/s 

w spanwise turbulent velocity, m/s 

w' spanwise turbulent velocity fluctuating about x — z average, m/s 

w s velocity of Stokes layer, m/s 

W wall velocity, m/s 

W m maximum wall velocity, m/s 

Wm^min minimum value of W m for drag reduction, m/s 

W t h threshold spanwise velocity used in the computation of S, m/s 

y vertical direction, m 

5 boundary layer thickness, i.e. wall-normal distance at which 
the mean streamwise velocity equals 0.99f/oo for boundary layer flows, m 

5 S Stokes layer thickness (8 a = \J , m 

6 momentum thickness for free-stream boundary layers, m 
v kinematic viscosity of the fluid, m 2 /s 

p density of the fluid, kg/m 3 

t x time- and/or space-averaged streamwise wall-shear stress, kg m/s 2 

t z space-averaged spanwise wall-shear stress, kg m/s 2 

+ indicates quantities scaled by inner variables, i.e. u T and v 
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Fig. 1. Schematic of the physical domain. 
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Fig. 2. R vw = dv'w' /dy + as a function of y + at different oscillation phases from 
the start-up of the wall motion (W+ = 18, T+ = 125). 
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Fig. 3. Drag reduction data from the present study (• and ■) and from QR (o) as 
function of the scaling parameter S. Refer to Section 4 for discussion on the data 
indicated by black squares. 
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Fig. 4. Drag reduction data as function of period of oscillation at fixed D m , for 
-D+ = 100 (•), 200 (A and A) and 300 (♦ and o). Solid lines represent the corre- 
sponding estimates based on (3) and (4). Refer to Section 4 for discussion on the 
data indicated by the open symbols. 
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Fig. 5. Optimum period of oscillation T^t t D as function of D^. The solid line rep- 
resents equation (5), while dots indicate the values of T^ ptD determined from the 
DNS data. 
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Fig. 6. Contours of the amount of drag reduction as function of W+, T + (top) 
and .D+, T + (bottom). Darker colours indicate higher drag reduction and contour 
increments are from zero (thick lines) by 5%. 
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Fig. 7. Contours of the percent net energy saving P ne t as function of W+, T + (top) 
and Z?+, T + (bottom). Darker colours indicate higher net energy savings, contour 
increments are by 1% and dashed lines are for negative values (only values for 
Pnet > —4% are shown). The thicker curves denote a null net energy saving. 
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9. Drag reduction data as function of S fro m three numerical stu d ies a t 
= 200. Result s from iBaron and Quadrio) (|l996h (•) and Ijl-Choi et aD (|2002h 



Fig. 

Re T 

(channel flow) (A). lQuadrio and Sibillal (l2000l ) (pipe flow) 
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Fig. 11. Drag reduction data as function of S from vario us numer i cal st udies. Results 
from I Jung et al.1 (|l992h : channel flow at Re T = 2 00 O. iNikitinl (l2000h: pipe flow at 



and r+ < 150 (o 



Re T = 147 for T+ > 150 ( 

layer flow) (o), and lMivake et al.l (|1997l ): channel flow at Re T 



Dha nak and Sil dl999l ^ (boundary 
150 (A). 
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Fig. 12. Drag reduction data as function of S fro m two experim e ntal s tud- 
ies on free-stream boundary layer flows. Resul ts from Ricco and Wu ( 20041 ) (o) 
at Re T = 257,451,633 for T+ < 150 and ISkandaiil (|l997h for T+ < 150 
(Re T = 375,463) (A) and T+ > 150 (Re T = 711) (a). 




Fig. 13. Dra g reduction data as function of S from three experimental studies. 
Results from Truiillo et al. ( 19971 ) (boundary layer flow) at Re T = 633 for T + < 150 
(biase d □ and corrected values o), T + > 150 (■). Results from IChoi and Graham 
(Il998h fpipe flow) at Re T = 650,962 for T+ < 150 (o) and T+ > 150 (•), and 
from lchoi et all (jl998h (boundary layer flow) at Re T = 549 for T+ < 150 (A) and 



T+ > 150 A 
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